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Abstract. We consider the one-dimensional Katz-Lebowitz-Spohn (KLS) model, 
which is a two-parameter generalization of the Totally Asymmetric Simple Exclusion 
Process (TASEP) with nearest neighbour interaction. Using a powerful mapping, the 
KLS model can be translated into a misanthrope process. In this model, for the 
repulsive case, it is possible to introduce second class particles, the number of which is 
conserved. We study the distance distribution of second class particles in this model 
numerically and find that for large distances it decreases as x~'^^'^. This agrees with 
a previous analytical result for the TASEP where the same asymptotic behaviour was 
found [T] . We also study the dynamical scaling function of the distance distribution 
and find that it is universal within this family of models. 



1. Introduction 

The study of stochastic particle systems with second class particles (SCP) [2j is 
twofold. On the one hand these systems provide relatively simple examples for particle 
systems with two conserved densities [3], since the dynamics of first class particles are 
independent of SCPs, while SCPs move stochastically on a landscape determined by 
the first class particles. On the other hand, the introduction of second class particles is 
often a powerful tool for studying fluctuations in the original model (with only first 
class particles) IH El [6], because SCPs in general follow the trajectories of density 
perturbations of first class particles. About the relationship between current fluctuations 
and fluctuations in the position of SCPs see e.g. [6J and [4j. Since SCPs are in general 
attracted by shocks which often appear in one-dimensional driven diffusive systems, 
SCPs can also be used for a good definition of the microscopic position of shocks 

PEimiHiEiiTnin]. 

In the present work we study the effective interaction between SCPs which is 
mediated by the common density landscape of first class particles. An important, and 
probably the most widely known example of a system where second class particles can 
be introduced is the one- dimensional Asymmetric Simple Exclusion Process (ASEP). 
The ASEP is a lattice model where each lattice site is either occupied by a particle or 
empty, double occupancy is forbidden. Particles can hop to the left and right nearest 
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neighbour sites with rates p and q {p ^ q) if they are empty [T^ [T3] . In this model 
second class particles can be introduced by the following dynamical rules [H |2l [T3] : 

10 — 7- 01 with rate p 

01 — i- 10 with rate q 

20 — )■ 02 with rate p 

02 ^ 20 with rate q ^ ' 
12 — )■ 21 with rate p 

21 — 7- 12 with rate g 

Here 1 and 2 indicates first and second class particles respectively, while is a vacancy. 
Notice that by not distinguishing between first and second class particles the dynamics 
in ([T]) reduces to that of the ASEP without SCPs. The same is the result if SCPs are 
not distinguished from vacancies. 

A breakthrough in the analytical study of the ASEP with SCPs came with the 
realization that the matrix product ansatz [HI [T5] can successfully be applied for finding 
the stationary state. In the seminal paper of Derrida at al. jl] the Totally Asymmetric 
Simple Exclusion Process (TASEP) with second class particles was considered (here 
g = 0). Among other important results about microscopic shock measures, it has been 
shown, that in the limit of zero density of second class particles, they form a weakly 
bound state. In particular, the distribution of the distance between two second class 
particles decays clS X 2 . 

We note, that related questions have been studied in [TB| [T7] for the Katz Lebowitz 
Spohn (KLS) model [18] with probe particles. The KLS model is an exclusion process 
with nearest neighbour interaction, which, for a special case, reduces to the ASEP, and 
in this special case the probe particles of [T6| [T7] are equivalent to the SCPs of the 
ASEP. Although the probe particles introduced in this model are in general not second 
class particles, the distance distribution between them was numerically found to decay 
as x~2 in [16]. 



3 

The purpose of this work is to check, whether the x~2 law is generic for second 
class particles of systems within the Kardar-Parisi-Zhang (KPZ) university class |19| . 
We consider a one-parameter family of misanthrope processes obtained by a mapping 
from the KLS model [20] and study the behaviour of second class particles in this process. 
In section |3] we introduce the model and define second class particles. In section H] the 
stationary properties (distance distribution of SCPs, density profile seen from a SCP) 
are studied, while section [5] is devoted to dynamical scaling properties. Results are 
summarized in El 



2. Second class particles in general 

Second class particles are most often studied in the context of the ASEP. Since other 
examples for SCPs are less well-known, a general introduction is given here. We consider 
systems where the occupation number rji of site i is integer, but is not necessarily 
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restricted to be between or 1. For a specific system rj^i^ ^ Vi ^ Vmax, where rj^ij^ and 
Vmax can be any integer, including ±oo. 

In general, second class particles (SCP) are defined through the basic coupling 
(described in detail e.g. in [6j) which works as follows. Two copies of the same system 
are considered simultaneously which are prepared in an initial state with occupation 
numbers that differ by one only at a single site. The time evolution of the two copies 
are then coupled in a way that marginally both of them satisfies the rules of the original 
model while at every time instance it is ensured that there is only a single site where the 
occupation number of the two copies differ. For models where this coupling is possible 
the mismatch site is considered to be occupied by one SCP. 

More than one SCP can be introduced in a similar way by considering two copies 
that differ at several sites with the condition that the following inequality holds for the 
occupation numbers 

Vf > Vf- (2) 

Here rjf denotes the occupation number of copy X at site i. The number of SCPs 
at site i is then defined to be rjf — rjf. Note that this definition is not restricted to 
exclusion processes, meaning that r/j, and therefore also rjf — rjf, can be any integer 
number. Due to the conservation law for the original particles, the number of SCPs 
Yli vf ~ vF ^1^° conserved. Systems, where property ([2]) at time zero implies that 
this property remains true for any t > are called attractive. In non-attractive systems 
the particle-interpretation of the occupation number differences is not justified, since in 
this case SCP - anti-SCP pairs can be created and annihilated (for the attractive case 
only annihilation is possible). 

Since marginally the two copies evolve according to the original dynamical rules, 
the presence of second class particles does not influence the dynamics of the original 
(first class) particles, so they can be considered as passive scalars. For details about 
the basic coupling see |6], where a wide family of attractive misanthrope processes are 
considered with several specific examples. 

3. Model 

3.1. One- dimensional KLS model 

The model considered here is the one-dimensional Katz-Lebowitz-Spohn (KLS) model 
|18] . which is a generalization of the Totally Asymmetric Simple Exclusion Process 
(TASEP). It is defined on a one-dimensional lattice of length L, where each site is either 
occupied by a particle or empty. During a continuous time evolution particles can hop 
to the right provided that the target site is empty. In contrast to the TASEP the rate 
of this jump depends also on the left and right nearest neighbour of the departure and 
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destination sites according to the following rules: 



0100 0010 with rate 1 

1101 1011 with rate 1 

0101-^0011 with rate 1-e 

1100 -> 1010 with rate 1 + e 



(3) 



Here 1 and indicates an occupied and an empty site. On top of the on-site exclusion, a 
positive e induces a repulsive interaction between particles, while negative epsilon leads 
to attraction. Note that here the sign of e is the opposite of what is usually used, i.e., 
e > corresponds to particle repulsion. 

For this system additional "probe" particles were introduced and the effective 
interaction between them mediated by the driven fluid (0 and 1) were studied in 
[I6| [T71 [2T| . These probes have been defined with the following dynamical rules: 



where 2 indicates a probe particle. In the special case of e = this model reduces to 
the TASEP with second class particles, where the second class particles are the probes. 
However, for any other e the probes do not behave as second class particles. In particular, 
the probes do influence the dynamics of Os and Is. 

It is relatively easy to see that the model defined in ([3]) is not attractive (in the 
sense described in section [2]) , the above described basic coupling would lead here to a 
system where SCP-anti-SCP pairs could be created and annihilated. However, as will 
be described in the next subsection, the model can be mapped to a misanthrope process, 
which is attractive (for e > 0), and where classical SCPs can be introduced (i.e., without 
SCP creation). 

3.2. Mapping to a misanthrope process 

The mapping, which has been used also in [20] consists of the following steps: 

• The vacancies of model ([3]) (i.e., sites with occupation number zero) are mapped 
to the lattice sites of the new model. 

• The length of the uninterrupted sequence of occupied sites to the right of a vacancy 
in the old model is considered to be the occupation number of the corresponding 
site in the new model. 

For an illustration see figure [1] As a result we obtain a model where the occupation 
number rji can be any nonnegative integer (r/min = 0, r/max = oo) with the following 
dynamical rules: 



12 — > 21 with rate 1 
20 02 with rate 1 



(4) 



rji, ^> r^j — 1, 1 with rate 1 + e 

Vi^ Vi+i ^ ^7* - 1, Vi+i + 1 with rate 1 

1, rji+i 0, r/j+i + 1 with rate 1 — e 

1,0 ^ 0, 1 with rate 1 



if r/i > 2 

if r^j > 2 and rji^i > 1 
if Vi+i > 1 



(5) 
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Figure 1. (a) Transition rates of the KLS model, (b) Transition rates of the 
corresponding misanthrope process. The rates depend on the occupation numbers 
of both the departure and the destination sites. 



This defines a misanthrope process where the hopping rate of a particle is a function 
of the occupation number of the departure and destination sites. It can be seen that 
this model, in contrast to the KLS model ([3]), is attractive for e > 0. (Note that 
here attractivity is meant in the sense described in section O the particles themselves 
interact repulsively for e > 0.) The main property that provides attractivity is that the 
transition rate r{rii,T]i+i), which is the rate from configuration r]i,rii^i to configuration 
rji — l, rji^i + 1 is a nondecreasing function of rji and a non- increasing function of r^j+i. For 
details see |6]. Notice that the particle- hole symmetry of the original exclusion process 
is lost by the mapping. 



4. Stationary properties 

In this section we consider two SCPs in the system with periodic boundaries and the 
distance between these two particles is studied. It is known that the dynamics of second 
class particles is closely related to density and current fluctuations in a system of the 
KPZ class [6], therefore results for SCPs reflect the nature of these fluctuations. In jl] it 
was shown analytically that in the ASEP the (non-cumulative) probability distribution 
of the distance between two second class particles in a stationary situation has an 
algebraic decay for large distances with a power 3/2: 

V{x) ~ x-t (6) 

(Here it is assumed that first the size L of the periodic system is taken to infinity, so 
1 < X <C L.) 

The question we study is whether remains true in a more general case. For 
this reason we study the misanthrope process introduced in section [3] with e > 0. Most 
studies of this kind assume half-filling for the exclusion particles which maps to density 
0" = 1 in the misanthrope process. We restrict our study for this special case. 

In order to check the above property we did numerical simulations in a periodic 
system with sites. Initially the system has been prepared in a state where the number 
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Figure 2. Shown is the measured distribution of the distance between two second 
class particles for various values of e in a periodic system of 3000 sites with a — 1. The 
result suggests that the exponent of the algebraic decay is universal (-3/2). Increasing 
the value of e the current-density relation becomes convex at p = 1. At this transition 
point the universality is not expected to hold 



of particles is with two extra SCPs (this corresponds to half filling in the exclusion 
model). Technically this means that we consider two copies in parallel: copy "A" holds 
N particles, while copy "B" is a clone of A with two extra particles, as described in 
section [2j Initially the particles are placed uniformly in the system, and after a sufficient 
relaxation time the distance distribution between the two SCPs is measured. Results 
are shown in figure [2j 

For most values of e it can be seen that the distance distribution follows ([6]), which 
suggests that this value is universal within the KPZ class. However, for a specific value 
of e (which appears to be near 0.9) a different behaviour is observed, which can be 
understood by looking at the stationary current-density relation. For systems of the 
KPZ class this relation is strictly nonlinear. The current-density relation j{p) for the 
KPZ model, and J(o") for the related misanthrope process is shown in figure [HI The 
mapping relates the two currents and densities as follows: 

where j and p are the current and the density in the KLS model while in the misanthrope 
process these are denoted by J and a correspondingly. An interesting property of this 
mapping is that the sign of the second derivative of j(p) always agrees with that of the 
second derivative of J{cf) at the corresponding value of a. 
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Figure 3. Shown is the current-density relation of the KLS model (a) and the 
corresponding misanthrope process (b). The functions j{p) and J(cr) are concave 
for e < 0.8. For half filling {p = 0.5, <t = 1) the second derivative of j{p) and J{<t) 
becomes zero at e = 0.8, and further increasing e it becomes positive. 



The analytical form of j{p) and therefore also of J(cr) is known (see [22]). It can 
be shown that the current- density relation is concave for e < 0.8, while for 0.8 < e < 1 
the second derivative becomes positive at p = 0.5 (a = 1). Near e = 0.8 and a = 1 
the second derivative vanishes which suggests a crossover to the Edwards- Wilkinson 
universality class, and this explains that the behaviour of V{n) is different at a special 
value of e (in the Edwards- Wilkinson class one expects no effective interaction between 
SCPs, which is consistent with the results of figure [2]). However, this value appears to 
be about 0.9 rather than 0.8, as suggested by the theory. It would be interesting to find 
a satisfactory explanation for this discrepancy. 

SCPs in general are attracted by high density gradients and the effective interaction 
between SCPs is mediated by the density fluctuations of the original particles. This 
can be clearly seen by measuring the density seen from a single SCP. Such density 
profiles are shown in figure HI Since the position of the SCP is highly correlated with 
density fluctuations, the profile measured from the SCP is in general far from being flat 
(although the stationary density profile, as measured by a standing observer is flat). 
The average speed of the SCP in a homogeneous environment with density cr is given 
by v{a) = J'{<j), which is also called the characteristic speed. At the point where the 
second derivative J"{cr) becomes zero the SCP becomes insensitive to (small) density 
fluctuations therefore one expects that the density profile seen from this particle becomes 
flat. This is what we see in figure |H Here the crossover point appears to be the same 
(ec ~ 0.9) as in figure [2] although the theory predicts ec = 0.8. 

Notice that even though for e = the KLS model reduces to the TASEP, the SCP 
of the corresponding misanthrope process is not equivalent to the SCP of the TASEP. 



Universal scaling for second class particles 



8 




Figure 4. The density profile as seen from a single SCP was measured in a periodic 
system of 3000 sites with a = 1 and several values of e (Ac = cr — 1). At a specific 
value ec the profile becomes flat which indicates a crossover to the Edwards- Wilkinson 
universality class. 



5. Dynamical properties 

It is interesting to see liow tlie distance between two SCPs clianges in time if they start 
at tlie same position. Similar studies for other sort of probe particles were done in 
|17) and [16]. We consider again the distance distribution between two SCPs but now 
as a function of time. In the simulation the system is prepared in a random initial 
configuration with a = 1 and two SCPs, just like in the previous section. After a 
sufficiently long relaxation time the two SCPs are monitored and the time is set to zero 
when the two particles occupy the same site (this initialization procedure agrees with 
the one used in [16j). The distance between the SCPs is measured and the distribution 
is calculated for sufficiently large number of independent realizations. 
Our scaling ansatz for V{x, t) for a specific e is 

V{x,t) = F[^)t-i-^, (8) 

where z is the dynamical critical exponent and F is a scaling function. The scaling 

3 

function F{y) for small arguments y is assumed to behave as F{y) ~ y~2 in accordance 
with the stationary behaviour. For the cumulative distribution function V{n,t) = 
YlT=n '^{^1 the above ansatz implies 

V{n,t)=G[^)t-^ (9) 
with another scaling function G, for which G{y) ^ y^^ ii y <^ 1. This scaling ansatz 
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10^3 lQ-2 iQ-1 iqO 

Figure 5. Shown is the scahng coUapse of the cumulative distribution function of the 
distance between two second class particles for various values of t and e. 



has been tested and confirmed in simulations. Furthermore, it is shown that not only 
the exponents but also the scaling function G is universal (up to non-universal scaling 
amplitudes a(e) and 6(e)). Results are shown in figure[5l where h{e)ti^V{n,t) is plotted 
against ^y/j with the dynamical critical exponent of the KPZ class z = 3/2. The 
amplitudes a(e) and 6(e) were fitted by hand for each value of e. We note that the 
scaling collapse is less convincing for larger values of e (e > 0.7) which might be due 
to the crossover expected at the point where the current-density relation turns from 
concave to convex. Close to this critical point the convergence is expected to be slow. 

6. Summary 

We considered a one-parameter family of asymmetric misanthrope processes, which are 
derived from the one-dimensional Katz-Lebowitz-Spohn (KLS) model [I8|. Although the 
standard KLS model is not attractive in the sense, that second class particles cannot be 
introduced with conservation, using a mapping to a misanthrope process, the resulting 
model (in a limited range of parameters) is attractive, thus second class particles can 
be defined using the basic coupling. 

We studied numerically the effective interaction between second class particles of 
this misanthrope process. It has been fond that the distance distribution between two 
second class particles has a power law tail with an exponent — |. This exponent has 
already been found to be exact in the Asymmetric Simple Exclusion Process [Ij. Note 
that although in a special case the KLS model reduces to the ASEP, the second class 
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particles of the misanthrope process are not equivalent to those of the exclusion process. 
It is interesting to compare this result with a related one for probe particles in the 
one- dimensional KLS model |16] . where similar distance distribution is found for probe 
particles, although these are not second class particles. 

In addition, the dynamical scaling function of the distance distribution has been 
obtained numerically. We found a good data collapse for several values of a model 
parameter, which suggests that this scaling function is universal within the KPZ 
universality class. 

For a special value of the model parameter e, the quadratic non-linearity in the 
(exactly known) current-density relation vanishes (at cr = 1), which suggests a crossover 
to another universality class (Edwards- Wilkinson). The crossover is observed in the 
simulation results, however, the location of this crossover shows a slight mismatch with 
the theory. This discrepancy could be due to finite size effects. 

The study of this paper could easily be extended to different densities of first class 
particles, and/or partially asymmetric dynamics. The behaviour for these more general 
cases is expected to be similar to what we found for the considered spacial case. Another, 
more interesting generalization would be to consider more second class particles and 
study the crossover between finite number and finite density of second class particles. 
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